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We present a theoretical study of Coulomb exchange interaction for electrons confined in a cylin¬ 
drical quantum wire and subject to a Rashba-type spin-orbit coupling with radial electric field. The 
effect of spin splitting on the single-particle band dispersioirs, the quasiparticle effective mass, and 
the system’s total exchange energy per particle are discussed. Exchange interaction generally sup¬ 
presses the quasiparticle effective mass in the lowest iranowire subband, and a finite spin splitting 
is found to significantly increase the magnitude of the quasiparticle-mass suppression (by upto 15% 
in the experimentally relevant parameter regime). In coirtrast, spin-orbit coupling causes a modest 
(1%-level) reduction of the magnitude of the exchange energy per particle. Our results shed new 
light on the interplay of spin-orbit coupling and Coulomb interaction in quantum-confined systems, 
including those that are expected to host exotic quasiparticle excitations. 
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I. INTRODUCTION 

The dimensionality of a many-particle system is a 
crucial determinator for how importantly interaction ef¬ 
fects can shape its physical properties. Generally, three- 
dimensional (3D) bulk conductors are less drastically af¬ 
fected by the Goulomb interaction between charge carri¬ 
ers than lower-dimensional, quantum-confined structures 
such as quasi-2D quantum wells and quasi-lD quantum 
(nano-)wiresii This is essentially due to phase-space re¬ 
strictions arising from free motion being only possible 
in fewer than three spatial directions. Furthermore, the 
exact structure of transverse bound-state wave functions 
shapes the density distribution of the confined charge car¬ 
riers and, thus, turns out to critically influence Goulom- 
bic effects in quantum wella^ and wires Here we explore 
how another aspect of quantum-confined states, namely 
their intrinsic spinor structure, modifies the effect of the 
Coulomb interaction in nanostructured systems. 

Most low-dimensional conductors are fabricated from 
semiconductor materials where the coupling between the 
spin degree of charge carriers and their orbital motion 
is often quite strong)^ As a result, quantum confinement 
can significantly affect spin-related propertiesi^ Such ef¬ 
fects are particularly pronounced for valence-band states 
(i.e., holes) because of their peculiar spin-3/2 charac¬ 
ter.—*^ In contrast, conduction-band electrons are spin- 
1/2 particles and generally subject to weaker spin-orbit 
couplings that are due to the bulk inversion asymme¬ 
try in the material’s crystallographic unit cell (Dressel- 
haus^ spin splitting) or the structural inversion asymme¬ 
try present in a nanostructured systems (Rashbcfc^ spin 
splitting). The multitude, and often counter-intuitive na¬ 
ture, of spin-orbit effects in nanostructures has become 
the focus of recent study, with developing an understand¬ 
ing of the interplay with Goulomb interactions being a 
key question to be addressed. Bulk-hole systeinsj^^— 
quantum-well-confined holeSfi^^— and 2D electron sys¬ 
tems subject to Rashba spin splittin g^^d^~ — have been 


considered. The comparatively few studies of Goulomb- 
interaction effects in spin-orbit-coupled quasi-lD sys- 
tems^“— have almost exclusively focused on effective 
Luttinger-liquid descriptions^! and, in particular, did not 
investigate the effect of Rashba spin splitting on the to¬ 
tal exchange energy and exchange-induced quasiparticle- 
effective-mass renormalization in quantum wires. 

In this article, we fill precisely this gap and investi¬ 
gate both the exchange energy and effective-mass renor¬ 
malization in quantum wires with a Rashba-type spin- 
orbit coupling. Previous work on the exchange energy 
of quantum wells revealed that spin-orbit coupling has 
the opposite effect on interactions in n-type and p-type 
systems: the exchange energy of a quasi-2D conduction- 
band electron system is slightly enhancedi^i^i due to 
spin-orbit coupling, whereas the exchange energy for 
quasi-2D holes is suppressed due to confinement-induced 
valence-band mixing!^ and the heavy-hole-type Rashba 
spin splitting!^ The different behavior of confined band 
electrons and holes warrants more systematic investiga¬ 
tion and, as we will see below, considering the quasi- 
ID case sheds new light on the different ramifications of 
spin-orbit coupling in interacting systems. Our investi¬ 
gation also reveals that the quasiparticle effective-mass 
is more strongly suppressed by the exchange interaction 
in nanowires with spin splitting. 

In addition, quantum wires with strong spin-orbit cou¬ 
pling are currently attracting great interest as possible 
hosts of exotic quasiparticle excitations such as Majo- 
rana3^ and fractional^! fermions. Clarifying the effect of 
interactions in such systems is necessary for a complete 
understanding of experiments aimed at verifying the ex¬ 
istence of the unusual quasi-particle excitations. 

The remainder of this article is organised as follows. 
We introduce our theoretical model of a Rashba-spin- 
split quantum wire in Sec. Ulland discuss pertinent prop¬ 
erties of the single-particle eigenstates. The formalism 
for calculating the exchange energy for this system is 
presented in Sec. lIIIi together with the results. Amongst 
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these is the ability to express functional dependencies of 
the exchange energy per particle in terms of a univer¬ 
sal scaling function, and the enhanced suppression of the 
density-of-states effective mass. Our findings are summa¬ 
rized, and related to the existing body of knowledge, in 
Sec.IVl Certain formal details are given in Appendices. 


II. THEORETICAL DESCRIPTION OF 
RASHBA-SPLIT NANOWIRE STATES 


In our study, we aim to develop a general under¬ 
standing of the effect of spin-orbit coupling on exchange- 
related many-particle corrections in quasi-ID nanowires. 
Hence, rather than attempting to describe the detailed 
electron density profile for a specific sample based on a 
self-consistent Poisson-Schrodinger calculation, we con¬ 
sider a model cylindrical quantum wire with radius R 
that is defined by a hard-wall potential where a constant 
radial electric field 6 = £r gives rise to a spin-orbit cou¬ 
pling of the Rashba type. In a real sample, such a radially 
symmetric field configuration could be generated, e.g., 
via biasing of an external gate that is wrapped around 
the wire surfacei^ The pragmatic assumption of a con¬ 
stant electric-field magnitude is justified in Appendix 
See especially Fig. [ 6 l For our situation of interest, the 
noninteracting-electron dynamics in the wire is described 
by the Hamiltonian H = -f U{r), where 


U{r) 


0 r < R 
QO r > R 


( 1 ) 


and is a Rashba-type^ single-electron Hamiltonian 


Here m* is the band mass of electrons in the semi¬ 
conductor material making up the nanowire, a is the 
material-dependent Rashba spin-orbit-coupling constant, 
and (T = {ax,cry,az)'^ denotes the vector of Pauli ma¬ 
trices. We will find the confined-electron states in 
the nanowire by superimposing solutions of the single¬ 
particle Schrodinger equation ip = Ep) io satisfy 
the cylindrical hard-wall boundary condition. 

The Hamiltonian (??) can be conveniently expressed 
in cylindrical coordinates (r, t/j, z) as 


i^(o) 


2m* 
+ ia£ 


1 d 

dj.2 j- Qj- 

1 d 

r o<f 


1 92 52 \ 


,(3) 


where a± = ((Ta, ± i (Ty)/2 are the spin-1/2 ladder oper¬ 
ators. The explicit form of ([3]) motivates a separation 
Ansatz for the eigenstates of : 


^ikz 


where 4>u,ki'r) is the radial spinor wave function, v = 
±1/2, ±3/2,... is an odd half-integer number, k denotes 
the wave number associated with the free electron mo¬ 
tion in the quantum wire, and L is the wire length. 
The resulting radial Schrodinger equation that deter¬ 
mines 4’u,k{f) can be written in dimensionless form as 
'Hu,kXuAq) = with 


m f I d\ lA _ rh ^ 

~ I Tr~^ H-I - 5 - OLGz - KGy-\-n 1, 

\oQ^ qoqJ Q 

(5) 

and the definitions m = u t \ Q = r/R^ k = kR, 
£ = E/Eq where Eg = h^/{2m*R^), a = 2Rm*a£/h^, 
and AxA = XyAf/R)/R- 

We employ the subband k • p metho d^^'^^ to find 
the cylindrical-nanowire eigenstates and single-particle 
subband-energy dispersions E^^J. Simultaneous invari¬ 
ance under time reversal _^ay = Hu^k) and 

spatial inversion Rv-k = 'Hv,k) imply that 

each subband is (at least) doubly degeneratei^ The first 
step is to find the eigenstates that are associated with 
the subband-edge energies A®. These states are then 
used as a basis set for expressing the eigenstates at 
general fc 7 ^ 0 ; with expansion coefficients determined 
from solving a matrix equation that is equivalent to the 
Schrodinger equation. 

The Hamiltonian of Eq. (??) is diagonal when k = 0, 


Riy,0 


H. 


H. 0 \ 

0 H_, y ’ 

/ 1_9 \ 

gdg) 


( 6 a) 


± 


(--hr 



, ( 6 b) 


hence the subband-edge states are also spin-projection 
eigenstates of cTj. with eigenvalue cr = ±1. We can there¬ 
fore write 


Xi^Ae) = I vA,n') ± I , (7a) 

n' — l 


with the subband-edge basis-state definitions 
I V, t, n') = AJA ^ (£•) (^ 0 ) ’ 

I V, 4 ., n') = ^ ^ , (7c) 


(e("b) 

and the functions {q) being solutions of the radial- 

confinement problem defined by the Hamiltonian + 
U{qR)/Eq with corresponding dimensionless eigenener- 
gies . We number the subband-edge states for fixed 
V and cr in ascending order of energy, that is > eA ^ 
when n' > n". Time-reversal symmetry mandates the 
Kramers degeneracy = AJ-a- See Appendix iBl for 
more mathematical details. 

The full single-electron subband dispersions can be 
found from solving the eigenvalue problem 
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FIG. 1. Electronic structure of noninteracting electrons in nanowires with spin splitting induced by a radial electric field. The 
solid curves in panel (a) [(b)] show the single-particle energy dispersions of the lowest two subbands obtained for a value of a 
corresponding to a recent experimental realization using InGaAs [InSb] as the wire material. To illustrate the effect of spin 
splitting, the corresponding dispersions for d = 0 are also plotted as dashed curves. Vertical lines are used to indicate the 
range of wave numbers for which only the lowest subband is occupied. Panel (c) illustrates more quantitatively the effect of 
spin-orbit coupling on the lowest nanowire-subband dispersions. In the upper (lower) panel, the ratio of the density-of-states 
effective masses (subband energies where SE = [E''°\a) — i5Q°^(a)]/[E^°^(0) — Eq°^( 0)]) for finite and for zero d are plotted as 
a function of wave number. 
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Ari'i) 


= £i^{k) 




( \ 

JU) 


In'i) 

V ; y 


(8a) 


with matrix elements 


(e) 




j(nn') ^ 27r / dp e 

Jo 

For the purpose of this study, we only need to obtain the 
dispersion of the lowest nanowire subband. We hnd that, 
for realistic values of a (see for instance the examples 
below), truncation of the eigenvalue problem (??) to the 
subspace spanned by the states {| 1/2, f, 1), | 1/2,4-, 1)} 
and its time-reversed counterpart yields sufficiently ac¬ 
curate results. Hence, in the following, we will use the 
wave functions 

^ikz 

ipi{r,ip,z) = -y^(- ismr]kR\l/2,t,l) 


^/LR 


^ikz 


e*‘^cos77fci^ 11/2,4,, 1)), (9a) 


'ip2{r,‘P,z) = - isinrjkR \ -1/2,4., 1) 

-I-e“*‘^cos77fc_R|-l/2,t, 1)) (9b) 
to describe lowest-subband states with the dispersion 


= Ei°} = Eo 


^Ik 


^2k 


(kRy 




^1/2,+ ^ 1 / 2 ,- 




( 1 ) 

1/2,-t 


-( 1 ) 
' 1 / 2 ,- 


+ (2dfci?/j)2^y 


.( 10 ) 


The coefficients entering Eqs. m are 


■ (8b) sm?7«,= 


1 

71 


cosr)^ = 


Ti) 

'1/2,-t 


-( 1 ) 

' 1 / 2 , 


v'(T/G-d‘,V)' + (2-A!;/)V 


(11a) 


1 

71 


1 - 


Ti) 
'1/2,-I- 


-( 1 ) 

' 1 / 2 ,- 


'W _M) 

^1/2,-l 1/2,- 


.-) +( 
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Figure [T] illustrates the noniiiteracting-electron band 
structure of nanowires using parameters relevant to re¬ 
cent experimental realizations^i in Ref. [3^ (InGaAs 
material with conduction-band effective mass m* = 
0.037 mg, where mg is the electron mass in vacuum, 
R = 300nm, and aE = 10“^^ eVm), and Ref. [s^ (InSb, 
TO* = 0.013TOg, R = 50nm, a£ = 10“^°eVm). Within 
our model, the relevant quantity determining the effect 
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finite k. Table |T] summarizes properties of the three low¬ 
est doubly degenerate subband edges in the two material 
systems. Note the rather large energy splitting of the 
(doubly degenerate) next-to-lowest subbbands due to the 
spin-orbit coupling. Without spin-orbit coupling (d = 0) 
the band edge energy of the subbands n = 3,..., 6 is 
^ 14.68. 

III. EFFECT OF SPIN-ORBIT COUPLING ON 
THE COULOMB-EXCHANGE ENERGY 


FIG. 2. The magnitude of the expectation value for spin 
projection parallel to the wire axis, Kuz)! = (ctz)! = —{pz)i = 
2 sin^ rikR — 1, for states from the lowest doubly degenerate 
(n = 1 and 2) subband. 


of spin-orbit coupling is a, which is equal to 1.82 and 
1.06 for the InGaAs and InSb nanowires, respectively. 
For comparison, we show also the result for d = 0. As 
the lowest subband-edge states have quantum numbers 
{v = 1/2, t} and {v = —1/2,4.}, respectively, their en¬ 
ergy is independent of d, and the spin-orbit coupling only 
affects the dispersion at finite k. In panel (c), the up¬ 
per plot shows the ratio of the single-particle density-of- 
states effective mass, 


mg (a, k) 


h?k 

dEf^/dk 


( 12 ) 


of the lowest subband with and without spin-orbit cou¬ 
pling for the two values of d. The lower panel in panel (c) 
illustrates the relative change in energy for the lowest- 
subband states due to the spin-orbit coupling. As can 
be seen from the plot, the renormalization of the single¬ 
particle effective mass due to spin-orbit coupling can 
amount to up to 30% (for d = 1.82) and also depends 
appreciably on the value of the wave vector. In Fig. [21 
we show the magnitude of the expectation value of the 
spin projection along the wire axis for the lowest sub¬ 
band as a function of the wave number k. It decreases 
with increasing k, as the states given in Eqs. (O together 
with m become superpositions of f and j. states for 


TABLE I. Properties of the three lowest doubly degenerate 
nanowire subband edges obtained for parameters applicable 
to recent experimental realizations. 


subband 
index n 

F/oVFo for 
a = 1.82 

for 

a = 1.06 

subband-edge 
(basis) state 

1 

5.783 

5.783 

+ ^,t, l) 

2 

5.783 

5.783 

|-il,l) 

3 

10.87 

12.47 


4 

10.87 

12.47 

1-14,1) 

5 

18.35 

16.85 


6 

18.35 

16.85 

+1,4, i) 


The Coulomb exchange interaction between electrons 
renormalizes the quasiparticle dispersion of nanowire 
subbands, which is then given byi 

(13) 

in terms of the non-interacting subband energy dispersion 
^nk obtained in the previous Section and the exchange 
(Fock) self-energy 

sS’ = - E/ ’•■’(■B"'*') ■ (w) 

n' 

Here n-p{E) denotes the Fermi-Dirac distribution func¬ 
tion, and V))/)” ^ is the matrix element of Coulomb inter¬ 
action between nanowire-electron states given by 


zr(nn') 

''kk' 


= C / d^r'L / 

J J J~L/2 

X Cl'fc'(l-_L)Cnfe(r_L) 


^i{k' —k)z 

dz — - 

\Jz'^ + |r_L - 

(15) 


where C = e^/(47reo£r) is the Coulomb-interaction 
strength, rj_ = (r, (p) denotes the position vector in 
the coordinates perpendicular to the wire axis, and 
Cnfe(i‘_L) = (t>u,k{r) is the transverse spinor 

part of the wavefunction in Eq. (??). In the follow¬ 
ing, we consider the zero-temperature limit and thus re¬ 
place the Fermi-Dirac distribution function by n-p{E) = 
&{Ey — E), with Q{E) being the Heaviside step func¬ 
tion and Ep denoting the Fermi energy. The condition 
Enk = Ef defines the Fermi wave vectors kpn for oc¬ 
cupied nanowire subbands. We now focus on the low- 
density situation where only states in the lowest doubly 
degenerate subband are occupied up to the Fermi wave 
vector /cp = kpi = kF 2 - For this situation, we can write 

Aintra(d, Kp , /t) 4- Ainj,ei.(d, Rp , ^) , 

(16a) 

where Ajntra (Ajnter) includes contributions arising from 
the exchange interaction between particles from the same 
band (from different bands). In the limit L —>• oo, we 
obtain the explicit expressions 
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sin^ sin^ r]^, 




fe'"' ) 
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■) 


+ cos^ ?7k cos^ 77k' 


J' 


■(4y2,-: 


(P) -Fir-'(770 


2 Ce(i) -) 2 

->.^1/2,-7 / /X 


/ "I 2 (rO'-') ) 2 fpti) -1 2 fgCi) -) 2\ 

+ sinaK cosOk sinaK' costtK'^ {g) [g') + Jq {g') T{g) j 


, (16b) 


^277 


Ainter(a,ACF,^) = / df<i' SlU^ {t}^ - p / dp'p' / d(^ COS 


/O ^0 


X {\k- K'\y/g‘^ + gi'2 - 200' cosi^) Jp (e) (g) Jp (s’') (&') , (16c) 


where Kq is the modified Bessel function of the second 
kindi^ For the numerical evaluation of the intra-band 
contribution (ll6bF we employ a modified quadrature 
method;^ described in greater detail in Appendix ICl to 
deal with the logarithmic singularity encountered when 
the argument of Kq(-) approaches zero. 

The exchange-renormalized density-of-states (quasi¬ 
particle) effective mass for the lowest subband can be 
calculated from 


mint{a,k) 


h^k 

dk ' dk 


(17) 


In Fig. [3l we compare the suppression of the quasipar¬ 
ticle effective mass due to the exchange interaction in a 
nanowire with finite spin-orbit coupling with that of an 
identical nanowire having zero spin-orbit coupling. As 
can be seen, the presence of spin-orbit coupling further 
suppresses the exchange-related quasiparticle-inass by 



FIG. 3. The ratio of the exchange-renormalized effective 
quasiparticle mass mint(Q;, k) to the bare single-particle ef¬ 
fective mass mo (a, k) is plotted as a function of wave vector 
k for states from the lowest nanowire subband assuming a 
material with dielectric constant Et = 12.9, kpR = 2.5, and 
spin-orbit-coupling strength d = 1.82 (fi = 0) as the solid 
(dashed) curve. For comparison, the dotted curve shows the 
result obtained under the assumption that the electric field 
strength varies linearly with the radial coordinate near the 
wire’s center, as described in Appendix lAl 


10-15% for the parameters used in our calculation. Note 
also the strong wave-vector dependence of the exchange- 
renormalized quasiparticle effective mass. 

The total exchange energy per particle for the 
nanowire-electron system is given by^ 


Ex 

w 


J_ ^ (X) 

2p^J 2 tt 


npiEnk) , 


(18) 


where p = N/L is the quasi-lD electron density. Again 
we focus on the low-density situation where only states 
in the lowest doubly degenerate subband are occupied up 
to the Fermi wave vector. For this situation, we can write 

Eyi C 

— 277 [■^intra(o^; T^f) T Ainter(d, T^f)] 7 (^^) 

where Aintra(a, kf) = f^^^dn Aintra(ci, kf, «:), and 
the analogous expression applies for Ainter- Figure 0] il¬ 
lustrates the functional dependences and relative magni¬ 
tudes of Aintra and Ainter- As Can be seen, the intra-band 
contribution is generally dominant and weakly dependent 
on a values considered here. In contrast, the inter-band 
contribution changes significantly as a function of a. 



FIG. 4. Scaling functions Aintra and Ainter associated with 
the intra-band and inter-band contributions to the exchange 
energy per particle in cylindrical nanowires with spin-orbit 
coupling (note the scale of 10“^ for the inter-band contribu¬ 
tion). Dashed (solid) curves corresponds to a = 1.06 (1.82). 
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FIG. 5. Relative change Ax in the magnitude of the ex¬ 
change energy resulting from a finite Rashba-type spin-orbit 
coupling quantified by parameter d, as defined in Eq. (??). 
Note the scale factor of 10~^ for the abscissa. The dashed 
(solid) curve shows the result obtained for a = 1.06 (1.82), 
which corresponds to a recent experimental realization using 
InSb (InGaAs) as the wire material. 


For quantum wires without spin splitting, i.e., in the 
case (5 = 0, the exchange energy per particle was found 
to obey a universal scaling formi^ i^^d° Our expression for 
Ex/N given in Eq. (??) generalizes these previous results 
to the case where spin-orbit coupling is finite. The change 
in magnitude of the exchange energy arising from finite 
a can be quantified through the relative difference 


_ Exja yf 0) 

Ex{a = 0) 


( 20 ) 


which is visualized in Fig. [5l For the values of a that 
correspond to recent experimental realizations using In- 
GaAs^ and InSb,— the associated change amounts to 
a suppression of the exchange-energy magnitude which 
can be up to 1.6%. This behavior is markedly different 
from the case of a 2D electron system where Rashba spin 
splitting has been shown^S to result in an increase of the 
exchange energy that is roughly one order of magnitude 
smaller. Thus the Rashba-type spin-orbit coupling due 
to a radial electric field in a cylindrical nanowire system 
is more similar to a 2D hole system where the interplay 
between quantum confinement and spin-orbit effects also 
results in a suppression of the exchange energy 


experimentally studied nanowire systems (one based on 
InGaAs and the other on InSb) as input, we have de¬ 
termined the influence of the spin-orbit strength on the 
lowest quasi-lD subband’s energy dispersion and on the 
spin projection of its corresponding eigenstates parallel 
to the wire axis, finding both quantities to be affected 
by tens of percent due to the presence of spin-orbit cou¬ 
pling. In particular, the density-of-states effective mass 
of the noninteracting system turns out to be increased by 
20-25% for parameters applicable to the InSb nanowires. 

With single-particle states in hand, we calculated the 
quasiparticle effective mass for the lowest subband and 
found its exchange-related suppression to be significantly 
larger in magnitude (by 10-15% for parameters used in 
our calculations) when spin-orbit coupling is finite. In 
contrast, the magnitude of the exchange energy per par¬ 
ticle is marginally reduced (by upto 1.6%) by spin-orbit 
coupling effects. Thus we find that any meaningful dis¬ 
cussion of the interplay between spin-orbit coupling and 
exchange interactions in quantum wires needs to be care¬ 
fully focused on specific physical quantities, as their rele¬ 
vant parametric dependences can be quite different, both 
qualitatively and quantitatively. Furthermore, often the 
relevance of interaction effects in an electron system is 
quantified in relative terms by a parameter rg that is re¬ 
lated to the ratio of contributions to the total energy 
arising from interactions and the single-electron disper¬ 
sion, respectivelyIn the present context, spin splitting 
causes an increase in the single-particle effective mass of 
quasi-ID electrons simultaneously with the suppression 
of the exchange energy. As the relative change in the 
increase in noninteracting system’s effective mass is an 
order of magnitude larger than the relative decrease of 
the exchange energy, the relative importance of interac¬ 
tions as measured by Tg turns out to be enhanced by 
spin-orbit coupling.— 

While we have focused on a specific configuration of 
confinement and spin-orbit coupling, our general results 
and overall conclusions can be expected to apply also to 
other spin-orbit-cc)upled nanowire systems, e.g., the one 
considered in Ref. [4^ 
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IV. CONCLUSIONS S. Chesi. 


We have studied theoretically the electronic properties 
of the quasi-lD electron system realized in a cylindrical 
quantum wire subject to a radially symmetric Rashba- 
type spin-orbit coupling. We determined the single¬ 
particle states for a hard-wall confinement using subband 
k • p theory. Focusing on the situation where only the 
lowest quasi-lD subband is occupied, we observed that 
the corresponding energy dispersion can be very accu¬ 
rately (to within 0.5% error) calculated from an effective 
2x2 Hamiltonian. Taking the material parameters of two 


Appendix A: Radial electric-field profile 

A proper self-consistent treatment of electrostatic ef¬ 
fects generally requires the application of an iterative 
Schrodinger-Poisson solver method that is specifically 
adapted to the sample layout. An added complication 
arises from the intricate way how the Rashba spin-orbit 
coupling strength needs to be determined from expec¬ 
tation values of the electric held taken in a multi-band 
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bound statei^ As we intend to focus on the broad im¬ 
plications of spin-orbit coupling in confined systems, we 
decided to make an assumption about the radial profile of 
the electric field entering in the spin-orbit term that en¬ 
ables us to obtain rather general physical insights. Here 

I 


2Trr L£(r) 



with the single-particle wave functions ip, z) 

given in Eq. Straightforward calculation yields 

£{r) = So [SkpR Voir/R) + Cu^r Vi(r/R)\, where So = 
—Ne/{2'KRL eo£r) is an overall scale containing the num¬ 
ber of particles N, = j: dn' sin^ 77 ^/ and = 
■i dn' cos^ ? 7 k' are weightings of the mixed bound-state 
contributions for the lowest nanowire subband, and 


^o(i)(£’) 



Q Jo 


) 



(A2) 


are the radial density profiles associated with the relevant 
bound states. 

The calculated full electric-field profile is shown in 
Fig. ini Our results from the main paper suggest that 
generally Sk^R ~ 1, Ck^R ~ 0; hence S{r) should be es¬ 
sentially determined by the Voir/R) contribution. This 
is indeed observed in the numerical evaluation. Also, 
as expected from the shape of the density profile associ¬ 
ated with the m = 0 bound-state wave function (cf. Ap¬ 
pendix [B]), the leading behaviour at r/i? <C 1 is linear. 



r/R 


FIG. 6. The radial profile of the electric-field magnitude S{r) 
associated with our calculated nanowire states is plotted as 
the solid red curve. The asymptotically linear behaviour for 
small r/R is captured by the blue dashed line. The horizontal 
(cyan) band indicates the range of field values that lie within 
10% of the mean value. The weak variation of S{r) for r > 
0.3 R motivates our assumption of a constant field magnitude 
for the radially symmetric Rashba term in Eq. (??). The scale 
of the electric field is eSo = 1.6kf /ieVnm“^ for the scenario 
based on InGaAs while it is eSo = 41kf /reVnm“^ for InSb. 


we show the basic consistency of this assumption with 
the electrostatics of the bound-state configuration for our 
system. 

Application of Gauss’s law using the cylindrical sym¬ 
metry of the nanowire geometry yields the relation 


dip / dr' r' 
Jo 


b/jir',P,z)f '4}jir',p,z) 


(Al) 


However, over most of the wire’s cross-section, the field 
profile is quite well approximated by a constant, which 
supports our pragmatic assumption. It is also observed 
from direct calculation that SkpR and Ck^R are almost 
constant in the relevant range k-pR < 2.5 where only the 
lowest nanowire subband is occupied. 


In order to confirm that the omission of the linear elec¬ 
tric field dependence for r/i? < 0.3 will not alter our con¬ 
clusions, we consider an electric field which is modelled 
by a linear dependence on the radial coordinate upto 
g < go and a constant for g > go- In terms of the di¬ 
mensionless Hamiltonian description in Eq. (??), this im¬ 


plies the replacement a ^ a 


j^Qigo - g) + Qig - go) 


Proceeding as in the case of a constant electric field, we 
find for the wave functions in the region with g < go 
the Bessel-function solutions + rna/go) ■ For 

the region g > go, we obtain wave functions which are 
a superposition of modified Laguerre functions (see Ap¬ 
pendix [B|) and confluent hypergeometric functions of the 
second kind. Applying the standard matching conditions 
at g = go to ensure continuity of the wave functions and 
their products with the velocity-operator in the trans¬ 
verse direction determines the unknown coefficients. In 
this context, it should be noted that the lowest state is in¬ 
dependent of the electric field. Considering the scenario 
with a = 1.82 as an example and taking into account 
the hard-wall boundary condition, we find only small 
changes for the band-edge energies E^^^-^^^/Eo = 11.07 

and = 18.22 when go = 0.3 (cf. Table [T|. In 

Fig. [7] we show the real part (dashed curve) and imagi¬ 
nary part (dotted curve) of the subband edge wave func¬ 
tion for spin up of the second-excited state | — 1 / 2 , f, 1 ) 
and compare this with the corresponding wave func¬ 
tion obtained under the assumption of a constant radial 
electric-field strength. We can therefore conclude that 
the linear electric field dependence for g < go changes 
the relevant wave functions used in our calculations only 
slightly. For k 0, we find that the matrix element 
~ —0.884, while it is ~ —0.916 with the as¬ 
sumption of a constant electric field, yielding only sub¬ 
percent changes for the dispersions and exchange-related 
quantities (see for instance Fig.[3|). 













FIG. 7. Real part (dashed curve) and imaginary part (dotted 
curve) of the subband edge wave function for spin up of the 
second-excited state obtained for an electric field that depends 
linearly on the wire’s radius upto g = 0.3 and is constant for 
Q > 0.3. For comparison we show the real wave function 
of the corresponding state (solid curve) obtained under the 
assumption of a throughout constant electric field. 


Appendix B: Solution of the radial-confinement 
problem 


The general solution of the differential equations 
present in the diagonal entries of Eq. (??) are power se¬ 
ries, given by, 

{n') 

Tt’^\Q) = g^{ao + aiQ + Y,an(r) (Bl) 

n—2 


which fulfill the relation {g) = {g) yield¬ 

ing the eigenstates. Disregarding the ill-behaved and un¬ 
physical part in the expansion at the origin, the coeffi¬ 
cients of the polynomials are determined by the recursion 
relation 


n{n ± 2m)an + maan-i + an -2 = 0 


(B2) 


with fli = —dm/(l ± 2m)ao-, where the upper (lower) 
sign applies to m > 0 (m < 0). The coeffi¬ 

cient uq is determined by the normalisation condition 

27r Jq dg g = 1. We note that the poly¬ 

nomial with coefficients given by Eq. (??) represents a 
modified Laguerre function that becomes the standard 

Bessel function Jo{\Jg) for d = 0 and/or m = 0. 

The band-edge energies, found by imposing hard 

wall boundary conditions on the radial wave function, i.e. 
for r = R, we require 


.Fi^'±>=l)=0. 


(B3) 


For not too large values of d, the lowest spin-f (i) 
subband-edge state has ly = 1/2 (—1/2) total angular 
momentum. However, as seen from Fig.jSl a level crossing 


FIG. 8. Energy eigenvalues of the lowest two doubly degen¬ 
erate quasi-lD subband edges, plotted as a function of the 
effective Rashba spin-orbit-coupling parameter a. The blue 
dashed curve is an approximation based on Eq. (??). 


occurs for d ~ 4.2, beyond which the new lowest spin-f 
(I) subband edge is a state with v = ±3/2 (—3/2). The 
variation of the band-edge energy , as a function 
of d can be approximated using standard perturbation 
theory, yielding 



- fo 

£q — a —j- , 

fo dg g Jlig^) 


(B4) 


where Eq ~ 14.68 is the band edge energy of the corre¬ 
sponding band for d = 0. 


Appendix C: Regularisation of the integrand for 
calculating the exchange energy 

In the calculation of the exchange energy we have to 
deal with integrals of the form 

1= jj dkdk'G{k,k')Ko {\k-k'Wr‘^ + r'‘^ - 2rr' cosLp^ , 

(Cl) 

with G(fc, k') being a smooth function of k and k'. A log¬ 
arithmic singularity occurs when the argument of Kof) 
vanishes. This happens when either the square root is 
zero, at f± = r'j_, or when k = k'. To regularise the inte¬ 
gral for the case where r_L = r'j_ , we add a small amount 
O'*" to the term under the square root. Then by decreas¬ 
ing the value of O'*", we perform a series of calculations 
until the result for the exchange energy doesn’t change 
within a certain tolerance. 

The situation for k = k' can be regularised analytically. 
To this end, we add to and subtract from Eq. (??) the 
term 


JJ dkdk'G{k,k)\n (^\k — k'\\/r^~+f^^'^^'2rR~cos^f'+fy^^ . 

(C2) 
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Adding this term to Eq. (??) cancels the logarithmic sin¬ 


J 


gularity. The fc'-integration of the subtracted term can 
be performed analytically and Eq. (??) becomes 


X = 


I I f^')Ko ^|fc — k'ly/r"^ + r'2 — 2rr' cost/?^ + G{k, k) In (^\k — + r'^ — 2rr' cosip^ 

kp + k ^ I _|_ j,i 2 _ 2rr' cos + 0+ 


-G{k,k) 


kp — k 


(C3) 


r 


The expression Eq. (??) is manifestly finite for k = k'. 
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